Nonequilibrium steady states of the isotropic classical magnet 
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We drive a d-dimensional Heisenberg magnet using a spatially anisotropic current of mobile par- 
ticles or heat. The continuum Langevin equation is analyzed using a dynamical renormalization 
group, stability analysis and numerical simulations. We discover a rich steady-state phase dia- 
gram, including a critical point in a new nonequilibrium universality class, and a spatiotemporally 
chaotic phase. The latter may be 'controlled' in a robust manner to target spatially periodic steady 
states with helical order. We discuss several physical realizations of this model and make definite 
predictions which could be tested in experimental or model lattice systems. 
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I. INTRODUCTION 

There is a growing awareness that spatio-temporal 
chaos may be a generic feature of driven, dissipative, 
spatially extended systems with reversible terms in their 
equations of motion 0, 0, OJ An understanding of this 
interplay between external drive, irreversible dissipation 
and reversible (inertial) dynamics may even be relevant 
to a deeper understanding of fluid turbulence described 
by the Navier Stokes equations. In light of this, a useful 
strategy is to construct simple spatially extended sys- 
tems whose dynamics reveals this interplay 0, Q and 
study its consequences in a systematic manner. In an ear- 
lier shorter communication we had successfully con- 
structed such a driven, spatially extended system which, 
while ordered in the absence of the drive, inevitably en- 
ters a spatiotemporally chaotic state when driven. In this 
longer paper, in addition to providing analytical and nu- 
merical details of our earlier work, we present some new 
results and discuss several microscopic realizations of the 
dynamics. We hope this study will encourage searches 
for possible realizations in lattice models and, more im- 
portant, experimental systems. 

The simplest examples of driven, spatially extended 
systems is the class of models called driven diffusive lat- 
tice gas models (DDLG) 0, . While these models give 
rise to a variety of interesting dynamics, none of them is 
chaotic. Recall that the DDLG order parameter is Ising- 
like; its dynamics is therefore purely dissipative and ex- 
ternally driven and has no reversible (inertial) terms in 
the equation of motion. The simplest way to incorporate 
reversibility is to elevate the Ising to Heisenberg spins; 
the driven dynamics of the 0(3) Heisenberg spins has a 
natural precessional dynamics [g. 

How do we drive the Heisenberg magnet spatio- 
temporally chaotic ? The dynamics of the classical 
Heisenberg magnet possesses several conserved quanti- 
ties (in fact there are infinitely many in dimension d = 1) 



. One possibility is to destroy some of these conserva- 
tion laws by driving it, not with an external field (which 
breaks the 0(3) symmetry explicitly) but by imposing a 
background steady current of heat or particles (or any 
other mobile species) which couples to the Heisenberg 
spins. If the dynamics of these mobile species is 'fast', 
then the imposed steady current may alter the effective 
dynamics of the isotropic magnet in a way so as to destroy 
some of the conservation laws of the un-driven system. 

In fact, as we saw in p], the classical 0(3) Heisenberg 
model in d = d± + 1 dimensions 8| in the presence of 
a uniform current in one spatial direction (|j), while re- 
taining isotropy in the order parameter space, should be 
governed by the equations of motion 

^ = (r«d? l +r 1 _V 2 ,)S-vS~-(S-S)S-\Sxd ll S 
ot \ 11 / 6 

+g ]] Sxdf l S + g x SxVlS + ff. (1) 

Eq. (0} is a natural generalization of the DDLG models 
to the case of a 3-component axial-vector order param- 
eter S and, as such, is an important step in the explo- 
ration of d yna mic universality classes | l0| far from equi- 
librium |3. The local molecular field in which spins 
precess in this driven state is dominated by the nonequi- 
librium precession term involving A, responsible for all 
the remarkable phenomena we predict, including a novel 
nonequilibrium critical point and, in a certain parameter 
range, a type of turbulence. Note that the noise rf in our 
equations is a non- conserving scalar noise even though 
the dynamics of the Heisenberg model is spin-conserving 
in the absence of driving. We shall comment below on 
the origin of this non-conservation as well as that in the 
deterministic terms. 

Here are our results : (i) Despite 0(3) invariance in the 
order-parameter space, the dynamics does not conserve 
magnetization; (ii) As a temperature-like parameter is 
lowered, the paramagnetic phase of the model approaches 
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a nonequilibrium critical point in a new dynamic univer- 
sality class controlled by the driving — we evaluate the 
dynamic z, roughening \ and anisotropy C exponents to 
leading order in e = 4 — d; (iii) Below this critical point, 
in mean-field theory without stochastic forcing, param- 
agnetism, ferromagnetism and helical order are all lin- 
early unstable; (iv) Numerical studies in space dimension 
d = 1,2, in the absence of stochastic forcing, show spa- 
tiotemporal chaos in this last regime. This chaos, when 
'controlled', is replaced by spatially periodic steady heli- 
cal states which are robust against noise. 

We also discuss a variety of explicit lattice and contin- 
uum realizations of this long length scale dynamics with 
the hope that this will stimulate a search for experimental 
systems, e.g., isotropic magnets carrying a steady particle 
or heat current, as well as model magnetized lattice-gas 
simulations, where the predictions of our model can be 
tested. Our work reinforces the idea that spatiotemporal 
chaos is a generic feature of driven, dissipative, spatially 
extended systems with nonlinear reactive terms. 

The generic occurrence of spatio-temporal chaos in this 
model is encouraging. It is instructive to compare this 
model with the dynamics of an incompressible fluid, em- 
bodied by the Navier-Stokes equations. Ignoring the fact 
that the while the total momentum of the fluid pv is 
conserved, the total spin in our model is not, the reader 
should not miss the formal similarity between the driven 
term S X duS and the nonlinear (v ■ V)v term in the 
Navier-Stokes, despite their disparate origins. The sim- 
ilarity of the convective nonlinearity of Navier-Stokes 
and the precession term in the dynamics of the classical 
Heisenberg model has already been remarked on @, 0] ; 
the form of our drive term (bilinear with one gradient) 
makes this resemblance closer. Our equations are how- 
ever more amenable to quantitative analysis, since, unlike 
Navier-Stokes, they are non-trivial even in d = 1. 

A section-wise breakup : In section II we provide 
a derivation of the equations of motion for the driven 
Heisenberg magnets (DHM) both from general sym- 
metry arguments and by analyzing explicit magnetized 
lattice-gas and continuum models. We next analize the 
steady state phase diagram of the continuum equations 
of motion (Section III), and show that the 'high drive- 
temperature' steady state of the model is paramagnetic 
(Sect. Ilia). There exists a nonequilibrium critical point 
whose critical exponents are determined using a dynami- 
cal renormalization group calculation (Sect. Illb). Lastly 
(Sect. IIIc) we study the 'low drive-temperature' phase 
of the system which exhibits spatio-temporal chaos. This 
chaotic phase may be 'controlled' giving rise to a steady 
state configuration with broken chiral symmetry (Sec- 
tion IV) . We end with a discussion of results and future 
directions (Section V). The details of the diagrammatic 
calculations are relegated to an appendix (Appendix A). 



II. DYNAMICS OF DRIVEN HEISENBERG 
MAGNETS 

Consider a system of Heisenberg (0(3)) spins situated 
either on a (hypercubic) lattice or a continuum in d- 
dimensions. Imagine a current of mobile species mov- 
ing along one spatial direction (say z), and interacting 
with the 0(3) spins. The spins interact with their neigh- 
boring spins via an exchange interaction (or any other 
short-range ferromagnetic interaction). If we treat the 
mobile species as being 'faster' (in a way which we will 
make precise) than the spins, then we may ask for the ef- 
fective dynamics of the Heisenberg spins themselves. We 
first derive the effective dynamics of the spins using gen- 
eral symmetry arguments and conservation laws. While 
such arguments are robust in themselves, they may not 
be realizable in a given physical setting. To counter this 
criticism we also construct explicit microscopic realiza- 
tions of the DHM. 



A. General Symmetry Arguments 

In this section we explain, how to construct, on general 
grounds of symmetry, the leading terms in the equations 
of motion of a uniaxially driven Heisenberg magnet. 

First let us remind ourselves of the equations of mo- 
tion for spins at thermal equilibrium at temperature T. 
The probability of spin configurations {Si} of a gen- 
eral nearest-neighbor Heisenberg chain with sites i is 
cx exp(— H/T), with an energy function 

H = — JjSj ■ 5,4-1, (2) 

i 

where Jj is the ferromagnetic exchange coupling between 
i and i + 1. A spin Si at i precesses as Si = Si x hi where 

-* OH -* 

hi = — = JiSi+i + Ji-iSi-i (3) 

dS z 

is the local molecular field. Replacing Jj — > J(x) and 

Si — > S(x) in the continuum limit, yields [l3j S(x) = 
J(x)S x dlS + (d.J/dx)S x d x S + .... For the physi- 
cally reasonable case where J varies periodically about 
a mean value Jo, this reduces for long wavelengths to 

S(x) — JqS x d^S, which is invariant under x — > —x, 
even if the H is not. If we impose a macroscopic dis- 
tinction between x and — x in the form of an (admittedly 
artificial) J varying linearly with x, say J — ax where a is 
a constant, we see that a term aS x d x S arises. However, 
the other term so generated, namely, J(x)S x 9^5, has 
a coefficient which grows with x, presenting problems in 
the limit of infinite system size. The dynamics in either 
of these cases conserves ^ Si since it commutes with H . 
Note that we have so far insisted on a local field arising 
from an energy function as in J5J|. 
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Let us now ask what the most general equation of mo- 
tion for S would be, if we relax the x — > — x symmetry 
and no longer insist on an energy function. This would 
clearly be the appropriate approach if the system were 
carrying a steady drift current of some mobile species - 
particles, vacancies, heat - in, say, the x direction, re- 
taining isotropy in spin space. Let us not specify at this 
stage how these additional degrees of freedom couple to 
the spins. It is sufficient to note that, irrespective of the 
nature of such couplings, the spins will be in a nonequi- 
librium state so that their dynamics will not follow from 
an energy function, and must be constructed anew. 

If we average over the degrees of freedom directly as- 
sociated with the current (the mobile species are 'faster' 
than the spins), their effect should be simply to modify 
the equations for the Si by allowing terms forbidden at 
thermal equilibrium. On general symmetry grounds the 
new terms will clearly involve an odd number of spatial 
derivatives. To lowest order in gradients, there are only 
two such terms: vd x S and XS x d x S. The first term, rep- 
resenting advection by a mean drift, may be eliminated 
by a Galilean transformation x — > x + vt, t — > t, leaving 
only the second term to reflect the drive. It is clear that 
the second term simply represents the continuum limit 
of asymmetric exchange, i.e., hi — J + Si+± + J-Si^i with 
A oc J + — J_ proportional to the driving rate. Such a 
term would be ruled out [8j at thermal equilibrium only 
because the dynamics had to be generated by J2J) and 10 • 

If we start with the usual Ma-Mazenko Q dynamics for 
S, with spin-conservation built in to both the systematic 
and noise terms, and then add in our novel precession (A) 
term, then standard perturbation theory for the noise 
and propagator renormalization yields, already at one- 
loop order, non-conserving terms of the form (S ■ S) n S 
in the limit of zero external wavenumber even if these are 
not put in at the outset. This is because the A term, while 
rotation-invariant in spin space, is not the divergence of 
a current 0,0] and, thus, docs not conserve total spin. 
A renormalization-group theory of the long-wavelength 
dynamics of such a driven system must allow from the 
start for such nonconserving terms. 

For a general dimension d = d± + 1, with anisotropic 
driving along one direction (||) only, the above argu- 
ments yield, to leading orders in a gradient expansion, 
the generalized Langevin equation (which possesses spa- 
tial 0(d— 1) symmetry), 

^ = (^r^ + r^l^S-vS-^iS-^S-XSxd^S 

+g l] Sxdf l S + g ± SxV 2 x S + ff, (4) 

where we have allowed explicitly for spatial anisotropy in 
all coefficients. The Gaussian, zero-mean nonconserving 
noise if satisfies 

(^(x, ^(x', t')) = 2B 6 af3 6 d (x - x')5(f - . (5) 

Note that the noise in our equations is a non- 
conserving scalar noise. It is vectorial only in internal 



(spin) space but is a scalar under spatial rotations. The 
correlation function of such a noise, at zero wavenumber, 
cannot detect spatial anisotropy. If the noise were con- 
serving, the covariance at wavevector q = (c/|| , qj_) would 
in general be proportional to gjj + aq]_, for some con- 
stant a 7^ 1. Such additive, conserving but anisotropic 
contributions to the noise doubtless exist here as well, 
but are irrelevant at small wavenumber relative to the 
nonconserving noise. 

Having established the form of the coarse-grained 
equations of motion on general symmetry grounds, we 
now offer explicit microscopic models within which the 
novel, nonequilibrium precession term could arise. 

B. Construction of Explicit Examples 

Recall that we are trying to create a driven state with 
directed spatial anisotropy (i.e., distinguishing, say, x 
from —x) while retaining isotropy in the spin degrees of 
freedom. Clearly, this cannot be achieved by means of a 
steady imposed spin current, since that would pick out 
a spin direction. We suggest here some approaches to 
achieve this. (I) Impose a steady current in some other 
species through the material and to show that this leads, 
via symmetry-allowed couplings between this additional 
degree of freedom and the spins, to the nonequilibrium 
precession term. The dynamics for the additional species 
is conserving] there is, strictly, speaking, no timescale on 
which the additional variables can be treated as fast and 
eliminated to yield an effective equation of motion for the 
spins alone. We must therefore employ a driven variant 
of model D, in the language of Ref. 10] . The effective 
spin dynamics we actually use in our paper should, how- 
ever, apply in the limit of an infinitely large diffusivity 
for the additional species or if processes that violate the 
conservation law for the additional species can be made 
to intervene. (II) Introduce the driving as a fluctuat- 
ing magnetic field, statistically isotropic in spin space, 
but with short-ranged correlations which distinguish x 
from —x. In both these cases (and, it seems clear, in 
general) it turns out that the exchange couplings have 
to be dynamical, not constant, in order to generate the 
nonequilibrium precession term. We now examine each 
of these cases. 

Consider a 1-dimensional lattice, each site i of which 
can be either vacant (ni(t) = 0) or occupied by one 
particle (ni(t) = 1) at time t. Each particle has an 
attached Heisenberg spin and may hop to the nearest 
neighbor at the right (left), if vacant, with probability 
p (q). The spin Si(t) at an occupied site i is the spin 
of the occupying particle. The local field at site i is 
hi(t) — Ji^ijSi-i(t) + Ji + i } iSi+\(t) where the exchange 
coupling Jij(t) determines, at time t, the field at j due 
to the spin at i. If we let depend on the configuration 
of the ni as Jij(t) = ni(t)nj (i)[Ji + (J2 — Ji)nj(t — 1)], 
with the Hi governed by an asymmetric exclusion process 
(ASEP) [a, [3 the driving nonlinearity in (0} is gener- 
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ated naturally. Explicitly, Ji,i±i(t) = J\ or J 2 accord- 
ing as rii(t — 1) = 1 or 0, since the exchange coupling 
between, say, i — 1 and i is operative at time t only if 
ni—i(t)rii(t) = 1, Assume for simplicity that the parti- 
cles can hop only to the right. Then a configuration f If 
at sites i — 1, i, andz + f at time t was either already 
present at time t — 1 or arose from Oil (by a right hop 
from i — 2). Thus, for f f f configurations, Ji-\s will be 
a weighted average of J\ and J2, while Ji+x,i will be J\. 
In the continuum limit, and averaging over the parti- 
cle dynamics, we will get the driving term in Q, with 
X rx p — q. 

The assumption of fast dynamics for the ASEP vari- 
able is justified if we allow evaporation-deposition as well, 
thus annulling the conservation law for particles. The 
derivation in our paper for the effective asymmetric ex- 
change felt by a spin at a given site remains unaltered 
by this non-conservation since, unlike the hopping, the 
evaporation-deposition affects in an unbiased manner the 
sites on either side of the spin in question. Of course, 
particle non-conservation induces spin non-conservation 
trivially in this case. 

More generally, let us construct a dynamical exchange 
coupling Jij(t) for a collection of spins in a host lattice 
which lacks invariance under x — ► — x. Imagine that the 
system is in a nonequilibrium environment, stationary, 
isotropic and translation-invariant in a statistical sense, 
in the form of a fluctuating magnetic field hi at every 
point i of the lattice. Assume that Jij(t) has a piece 
AJij(t) oc hi(t) ■ hj(t — t) where r is some timescale in- 
trinsic to the material. Since the material lacks x — > — x 
invariance, the average (AJy) cx (hi(t) ■ hj(t — r)) is 
in general nonsymmetric in ij if the sites i,j are sepa- 
rated in the x direction. Microscopically, this could hap- 
pen if there is some underlying dynamical process, for 
instance structural relaxation, determining the exchange 
couplings between sites i and j. Note that we do not need 
explicitly to drive anything through the system; asym- 
metry under x — > — x and a nonequilibrium noise should 
suffice to produce this ratchet-like effect. In addition, the 
precession of spins in this fluctuating field implies a term 
Si x hi which is not spin conserving. 

Curiously, a term with precisely the form of our 
nonequilibrium precession term appears in the equilib- 
rium dynamics of the staggered magnetization in the 
isotropic antiferromagnet, as presented in the work of 
[l3|. The consequences of such a term appear not to be 
very dramatic there, presumably because of the presence 
of other couplings to the ferromagnetic order parameter. 
In addition, equation J3J with no noise and with all co- 
efficients on the right-hand side except A set to zero is 
known in the literature as the Belavin-Polyakov equa- 
tion 0, 0| and has been widely studied for its soliton 
solutions. 

Moving to a completely different description, consider 
the dynamics of chiral self-propelled particles de- 
scribed by a polar vector n, suspended in a fluid and 
drifting along the z direction. A transient variation of 



n along z must lead |l9| to an overall rotation about 
ft. Shifting away the effect of the mean drift this clearly 
yields a dynamics of the form dffi ~ n x d z n + . . . which 
is once again our nonequilibrium precession term. The 
consequences of this for the dynamics of self-propelled 
chiral objects remain to be worked out. 



C. Fokker Planck Equation and absence of an FDT 

Since in the Langevin dynamics of this driven Heisen- 
berg magnet, both the nonconservative noise and deter- 
ministic terms arise from the external drive, there is no 
obvious relation (such as the fluctuation-dissipation the- 
orem (FDT)) between the variance of the nonconserva- 
tive noise and dissipation [2(| . This may easily be seen 
by constructing the Fokker-Planck equation [21j for the 
probability distribution of spins P({S},t) corresponding 
to the Langevin equation (J3J, 



dP 
~dt 



d 
dS a 



A a P — B 



dP 



dS a 



(6) 



where 



A = [r ll df l +rx^i)S-vS--(S-S)S-XSxd {l S 



gnS x e%S + g±S x \7\S 



(7) 



Since the stationary probability distribution of the steady 
state configurations need not be the equilibrium canon- 
ical distribution exp(— F[S]/T) at a temperature T 
(where F is the free-energy functional), there is no di- 
rect relation between the B and the other parameters 
that enter the Langevin equation. 



III. NONEQUILIBRIUM STEADY STATES 

Having discussed several microscopic realizations of the 
coarse-grained continuum dynamics Eq. |4"|. we will pro- 
ceed to establish a 'nonequilibrium phase diagram' of 
steady states obtained by analyzing the stationary solu- 
tions (d t S — 0) of Eq. (@J. Let us specify the parameters 
in our 'phase diagram'. 

In the equilibrium, isotropic limit, A = u = v = 0, 
ru = rx = r, the noise strength (which is conserved in the 
absence of the drive) vanishes at zero wavenumber, and 
Eq. has a critical point where the renormalized r goes 
to zero. In the driven state, the dynamics and noise are 
nonconserving; the critical point is at v = 0, which in gen- 
eral takes place on a curve in the temperature/driving- 
force plane. As the drive is taken to zero there should be 
a crossover from nonequilibrium to equilibrium critical 
behavior. Our primary interest is in the behavior at a 
given nonzero driving rate, for which it suffices to vary 
the temperature-like parameter v (which we shall refer 
to as 'drive-temperature') in J3J, keeping the rest fixed 
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(with ry, r±,u > 0). Note that the variance of the non- 
conserved noise B is another temperature-like parameter. 
Since the FDT is violated by the external drive, these two 
temperature-like variables v and B are not related to each 
other. The nonequilibrium phase diagram will therefore 
be parametrized by the two parameters v and B. 



A. Dynamics at High Drive- Temperatures (v > 0) 

At high drive-temperatures the steady state, obtained 
by setting dtS — is paramagnetic, characterized by 
(5(x, t)) — ((■ ■ •) denotes an average over the noise if) 
and correlators . The effec t of th e drive is to change the 
correlation lengths £j_ = ^r±/v + 0(X 2 /^/v) — 0(X 2 ) and 
£|| = v/^w + 0(X 2 /^/v) - 0{\ 2 ). Thus for v > (the 
paramagnetic phase) all correlations clearly decay on fi- 
nite lengthscales ~ 1/\A> an d time-scales ~ 1/v, and 
nonlinearities are irrelevant. We find that this paramag- 
netic state is stable under dynamical perturbations. This 
can be seen by writing S(x, t) = (5(x, t)) +u(x, i), where 
u is an arbitrary small perturbation. The time evolution 
of u(x, t) to linear order is given by 

du / \ 

— = [rxVi+r^-vju + ff, (8) 

which on Fourier transformation reads 

■u k (i) = u k (0)exp(-7 k t) 

+ dt'rk(t')exp(-y k (t-t')), (9) 
J o 

where 

7k = rj_k 2 ± + r\\k 2 + v . (10) 

This arbitrary perturbation always decays to zero 
when v > 0. The 'paramagnetic phase' is therefore lin- 
early stable. 



B. Dynamics in the Critical Phase (v — 0): New 
Critical Behavior 

A description of the nature of correlations on the crit- 
ical surface v = 0, requires the machinery of the Dy- 
namical Renormalization Group (DRG). While there are 
several clear expositions of this technique 22] devised for 
specific problems, we hope that some readers may bene- 
fit from a pedagogical treatment of DRG applied to our 
anisotropic driven Heisenberg dynamics. 

Let us first drop all nonlinear terms from the equa- 
tions of motion. In the critical region, defined by v = 0, 
the linear theory is massless resulting in divergent long 
wavelength fluctuations, as can be seen by explicitly cal- 
culating the correlation function C(x, t) = (5(x + x , t + 



t ) ■ S(x , t )) from Eq. when u = A = 0. The corre- 
lation function is easily recast as a scaling form, 

C( X ,t) = ^F ( r ^J^), (11) 

r ll V x \\ x>, \ r ±/ 

where the roughening exponent \ = 1 — d/2, the growth 
exponent z = 2, the anisotropy exponent ( = 1, and Fq 
is an analytic function of its arguments. 

What is the nature of these divergent fluctuations in 
the presence of the nonlinear terms? This can be ad- 
dressed via a standard implementation of the dynamical 
renormalization group (DRG) |22j based on a perturba- 
tive expansion in the small couplings A and u, about the 
linear theory. The perturbative corrections to the cor- 
relation function may be equivalently viewed as arising 
from modifications (renormalization) of the parameters 
ru, r±_ and B. Renormalizability guarantees that the cor- 
relation function C(x,t) will retain a scaling form as in 
Eq. Ullj l with modified exponents z, Q, \ an( i a new 
scaling function 

cm) = *;*f(±,|). ( i2) 

This implies that the critical region, defined by (renor- 
malized) v R — 0, still has divergent long wavelength fluc- 
tuations. In what follows we assume renormalizability, 
which we justify post facto to lowest order in perturba- 
tion. 

The perturbation expansion would make sense only if 
the couplings A and u are small. This is ensured by the 
DRG procedure which involves a systematic expansion 
about the upper critical dimension d c . As we shall see, 
the renormalized couplings A and u will flow to small 
values, of order e = d c — d. The upper critical dimen- 
sion may be obtained by a simple power-counting ar- 
gument. Rescale space, time and the order parameter 
as : X|| = bx'p x^ = b^x.' ± , t = b z t' and S = b x S', 
where b > 1 is an arbitrary parameter. We may rein- 
terpret the effect of such a rescaling as a change in the 
parameters ; thus the form of Eq. will remain un- 
changed if we change the parameters to their primed val- 
ues rf, = b z -' 2 r h r'^ = b z -^r ± , B' = ^- 2 x-C(d-i)-i 5) 

u' = b 4 - d u, A' = fcX+^A, g'n = b* +z - 2 g\\, and g' ± = 
b x+z ~ 2 g±. Demanding that the linear equation (in the 
absence of nonlinearities, u — A = g\\ = g± = 0) be scale 
invariant, automatically fixes x = 1 ~ d/2, z = 2, Q = 1, 
consistent with Eq. Qlip. 

Inserting the values of the exponents evaluated within 
the linear theory, we find that the couplings A and u 
are relevant, i.e., they grow under rescaling, for dimen- 
sions d < d c = 4, the upper critical dimension for these 
couplings. Thus nontrivial exponents are expected for 
d < d c = 4 in the presence of nonlinearities. The mean- 
field estimates in the previous paragraph are valid for 
d > 4. The other two couplings gu and g± have an upper 
critical dimension d 9 c = 2; they are therefore irrelevant 
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for d near 4, and will be set to zero in the analysis that 
follows. 

We now carry out the DRG calculation to compute the 
corrections to the new exponents and the scaling func- 
tion due to the nonlinear terms using a field theoretic 
approach [H HJ . The DRG procedure involves, 
as usual, two steps — 

(i) Start with a large-wavenumber cutoff A. Solve 
the equation of motion iteratively in an expansion in u 
and A Average over modes with wavenumber k in the 
shell (A/6, A), with b = et, to obtain "intermediate" 
renormalized parameters r^, r|, X 1 etc. in the equations 
of motion. This averaging is carried out to leading 
order in A and it, justified post facto by the fact that 
these couplings flow to small values, of order e, at the 
DRG fixed point in dimension d = 4 — e. (ii) Rescale 
space, time, and fields as in the previous paragraph to 
restore the original cutoff A. The parameters rj_ etc. 
will then acquire rescaling factors as above. Define 
r±(£) = fo z ~ 2< »r^, and similarly define ^-dependent 



versions of the other parameters in the equation. This 
may be recast as differential recursion equations in the 
various couplings. 

(i) Perturbative Calculation : The bare (unrenormalized) 
propagator Go(k, to) and the correlator Co(k, to) are de- 
fined from the linear theory as 

G (k,cj) = 2 ^2 — i — > 

7~|| «j| + r±kj_ + v — ltd 

C„(k,w) = (£(k,w)-£(-k,-w)) 



(r\\k?, + r±k\ + v) 2 +lu 2 ' v ' 

We next calculate the corrections to the correlation 
functions from the nonlinearities, perturbatively in the 
couplings u and A. On Fourier transforming Eq. Q we 
obtain 



iX f 

S a (k,uj) = G (k,to)r) a - yG (k,w) / 3qSv[q\\ - {k\\ - q\\)]e a psSp{ci,v)S 5 {k- q,w - v) 
-uF a p lS Go(k,uj) J dq 1 dq 2 diy 1 diy 2 Sp(c[ 1 , v 1 )S~ l {q 2l v 2 )Sp{k - q x - cfe, to - v x - v 2 ) , 



where the Fourier transform is defined as 

S(x,t) = J dkdu;S{k,uj)e- luJt + lk - x , (15) 

with the measure dq = d d q/(2ir) d , dv = dv/2n and 
the range of integration < \q\ < oo, — oo < v < oo. 
The coefficient of the cubic term is given by F a f3 7 s = 

(l/3)(8 a p8 7 8 + SatSpy + 5arf8ps)- 

Note that we have explicitly retained v in the above 
expressions; after we have calculated the renormalized 
correlators, we shall set the renormalized v R = 0. This 
defines the critical phase. To carry out the perturba- 



(14) 



r 



tive calculation effectively it is convenient to rewrite the 
recursion relation Eq. (|14|) in terms of the graphical rep- 
resentation displayed in Fig.^ The graphical representa- 
tion is standard — with — * denoting the bare propagator 
Go(k, to) and x denoting the noise r) a (k, to). The averag- 
ing over the noise is performed using 

(r ]a (k 1 ,uj 1 )r]f i (k 2 ,uj 2 )) = 2BS d (k 1 +k 2 )S{io 1 + Lo 2 )S af3 . 

(16) 

The renormalized propagator G R (k, to) can be ob- 
tained perturbatively from Fig. ^ The lowest order (1- 
loop) correction is given by, 



(<7|| - (h -q\\))(k\\ - (-<7||))G (k-q,w-i/) 



G R (k,to) =G (k,w)-Go(k,w) 2 j^- J dqdv 

xG (q,z/) + ((fc|| -g||) +(fc|| - g||))G (q,^)Go(k-q,w- v) 

+y J dqdvC {(\, i/)\ 



Go(k,w) + G (k,w) 2 £(k,w). 



(17) 



Note that the 1-loop renormalized propagator has terms to 0(A 2 ) and 0(u) in the couplings. The self energy 
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i ff. q ■ v ) 



T) ( k ,(0 ) 



r\ ( k , ra ) 



S„(q,v) 



S a (k,co ) G„(k,to ) G (k,(o) 

S y ( k-q , co-v) 11 



Tlyf k-q , CO -v) 



\(q,.v,) i,(q 2 .v 2 ) 



Sp(q rVl ) 



s y (q 2 .v 2 ) 

S 5 (k-q-q [O-v -V ) 



G (k,to) ll 6 (k-q-q 2 , co-v r v 2 ) 



(a) 



(b) 



^\[q H - (k„- q„)] 



j 

q 



uF a|3yS 



q,,q 2 



FIG. 1: Graphs depicting the recursion relation in S. A thick 
line with a cross at one end denotes S(k, ui) and a thin arrowed 
line denotes Go(k, lo). The cross x denotes the noise r/ a . 
(a) and (b) denote the vertices coming from the A— and u— 
couplings respectively. 



ten as 



[G^k,^]- 1 = r R k\+r R kl + v R - 



(18) 



to obtain the renormalized parameters v R , r R and r R . 
Thus in order to read out the corrections to the bare 
parameters, we need only compute the lo — > limit of E, 
and expand the 1-loop corrections to E in powers of the 
external momenta, retaining only the coefficients of k°, 
k 2 and k\ (Appendix A). 

The 1-loop corrections to the parameters are given by, 



-6/2 



4e 



7 ± 



r± 1 + 



5lTTV £ / 2 

"48e 



7T Ar ± - + 



5 K 

"6~ 



(19) 



(20) 



(21) 



E(k, lo) defined above contains all the 1-loop corrections 
coming from the nonlinear terms. The calculation of 
E(k, lo) at the critical point (v R — > 0) reduces to evaluat- 
ing integrals — these are singular in the infrared k — > 
limit for d < 4. This apparent divergence reflects the 
relevance of nonlinearities below 4 dimensions; indeed it 
is this kind of divergence that the renormalization group 
procedure is geared to handle. For now we note that the 
integrals turn out to be logarithmically divergent, both 
in the infrared and the ultraviolet k — > oo limit, as d — > 4; 
further the divergent pieces in both the limits turn out 
to be the same [23). This allows us to use a procedure 
known as dimensional regularization [25^ . to isolate the 
divergences as poles in 4 — d = e = 0. The ultravio- 
let divergences are controlled by a momentum cutoff A. 
Details of the calculation are presented in Appendix A. 
The 1-loop renormalized propagator G R may be writ- 



where the effective coupling constants r and k are defined 
as 



X B 



and 



uB 

r, 1/2 3/2 



(22) 



(23) 



As discussed the ultraviolet divergences have been regu- 
lated by an upper momentum cutoff A while the infrared 
divergences appear as poles in e (dimensional regulariza- 
tion) . 

To evaluate the corrections to the other parameters we 
need to calculate the renormalized correlation function 
C R (k, lo) = (S^k, lo)-S(— k, — lo)) which to 1-loop satisfies 



C R (k,Lj) =C (Klo 



4 



Sqdv({k\\ - q\\) - q\\){q\\ - {k\\ - g||))C (k - q, w - v)C {^v) 

xG Q (k,Lo)G Q (-k,-Lo) 
= 2G (k, Lo)G Q {-k, ~lo)B + G (k, u)G Q {-k, -lo)Z b . 



(24) 



The function E^ contains all the corrections coming 2B R G R (k, lo)G r (— k, — lo), 
from the nonlinear couplings and is calculated in Ap- 
pendix A. The renormalization of the noise may now / nTV -e/2 \ 

be easily determined via the definition C R (k, lo) = B =51 1H — I. (25) 
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We have also computed the lowest order corrections 
to the nonlinear couplings (vertex corrections) defined as 
— i\ R — — iA+2 T\, and — u R = — u+T u , where the vertex 
functions T\ and T u have been computed in Appendix A. 
To 1-loop we find 



X H = X 1- 



3tttv~^ 2 
128e 



and 



u — (' 



1/2 3/2 

*ll r -L 



ll7T 4 « 2 

12Be 



27^ 4 r 2 
8Be 



(26) 



(27) 



(ii) Recursion Relations : The relation between the renor- 
malized and bare couplings allows one to write down the 
scaling behavior of these couplings. First note that the 
parameter v appears in the propagator as a mass term; 
infrared divergences characterizing critical behavior arise 
only in the massless theory. Trading this mass scale for 
a scale of length, we may define an "observation scale" 
£ = bao ~ b/A = v~ x l 2 , where do is a microscopic cutoff 
length and b > 1 is a pure number. Thus, large k/v 1 / 2 
means small x/£ which defines the critical region. 

Of course if the limit b — > oo is taken straightaway, the 
renormalized couplings would diverge. The strategy is 
therefore to increase b gradually averaging the couplings 
over a small interval between b and b+db — this amounts 
to writing a differential equation for the evolution of the 
couplings as the scale parameter b changes. 

In the critical region, one may define dimensionlcss 
couplings using the above observation scale; for instance, 



(b) 



l (ba y 



As may easily be seen from Eq. 



(|19fl . this dimensionless coupling satisfies, 



f R 



(b)=r l{ (ba y 



1 



TTT(ba ) 
4e 



Applying the rescaling operator bd/db, we obtain, 
df R 



II 

db 



r\\(ba ) 



TTT(ba o y 



(28) 



(29) 



Similarly define a dimensionless parameter for r; f (b) = 
T R (ba y. Recalling that r is of 0{\ 2 ) (Eq. (JUJ), we 
may to lowest order in r replace and r in Eq. 1|29|) 
by their renormalized values. Finally expressing in terms 
of I = In b, we arrive at the 1-loop differential recursion 
relations (the ultraviolet cutoff A has been set to 1), 



d 



dt 
8rf 
~dl 
dB R 

d£ 
0f R 

d£ 
dk R 
~d£ 



f R (z- 



f R ) 



= ? R (z-2( + 



■ r )7T 



f R ) 



B 1 



f R 



(Ce 



48 

2 X -((d- 

64 h 
11 
24 



32 



27 
16 



Recall that the couplings A and u enter the per- 
turbation theory in the dimensionless combinations 

r = (1 / '2tt 3 )A 2 Bj 'yVjj r\ and k = {1/2tt 3 )uB/ 

The reader would have noted that we have set the 
irrelevant couplings g« = g± = for d = 4 — e at the 
critical point, v R = 0. 

(in) Fixed Points and RG flows : We expect the scaling 
behavior Eq. to hold at the critical point in the limit 
£ — > oo. Thus the parameters f R , etc., should be scale- 
independent for £ — > oo. This necessarily implies that 

-gj- = as £ — » oo and so on. In other words, the critical 
behavior at v R — is given by the fixed points of the 
recursion equations derived above, 



Of* 
~d£ 



dk* 
~d£ 



~d£ 



drr_d_X L _d_B 1 _ 
d£ ~ d£ ~ d£ {6 > 



The above equations yield four fixed points; by mak- 
ing small deviations from these fixed points along the 
remaining directions in parameter space we can deter- 
mine their stability (to linear order). The exponents z, 
£ and x mav be evaluated (to 0(e)) at the stable fixed 
points : 

(A) f*—k*=0 is the 'gaussian fixed point' and is sta- 
ble for d > 4 and unstable for d < 4. The expo- 
nents take their 'mean field' values z = 2, £ = 1, 
and x — 1 — d/2 at this fixed point. 

(B) f* — 0, k* — 24e/ll7r is an unstable fixed point 
both for d > 4 and d < 4. 



(C) f* = 646/357T , k* = 24e[3 + V12235] /385tt is an- 
othcr unstable fixed point both for d > 4 and d < 4. 



(D) f* = 646/(35tt), k* = 36e [1 + yT409]/3857r is a 
nontrivial 'driven fixed point' and is unstable for 
d > 4 but stable for d < 4. Exponents take non- 
trivial values to 0(e), z = 2-16e/35, C = l-2e/15 
(anisotropic), and x = 1 — d/2. Note that x does 
not change from its mean field value to this order, 
since the quartic vertex u plays no role at lowest 
order in e. 

Thus for e = 4 — d > 0, the nontrivial stable fixed point 
(D) is associated with the critical exponents z = 2 — 
16e/35, C = 1 — 2e/15 and x = 1 ~ d/2, to lowest order in 
e. These exponents clearly place this critical behavior in 
a new universality class and different from the anisotropic 
KPZ 26] In Fig. |3 we exhibit the fixed points and the 
RG flow diagram to 0(e) in a 2-dimensional subspace of 
the entire parameter space. 

We end this section with a brief remark. There has 
been some recent discussion in the literature 0, Q that 
for a particular class of driven 0(n) models, there exist 
no stable fixed point to given order in e. The authors 
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FIG. 2: Plot show fixed points and RG flows. For d > 4, A 
(Gaussian fixed point) is the only stable fixed point, while for 
d < 4 there is a nontrivial stable driven fixed point D. 



suggest, without proof, that the RG flows in these mod- 
els may be chaotic. We would like to stress that in our 
model, there exists a stable fixed point in all dimensions, 
and so the critical behavior here, atleast to O(e), is not 
chaotic. However we will see in the next section that 
the entire low drive-temperature phase, in the absence 
of noise, is chaotic. We will get back to a discussion of 
these issues towards the end. 



IV. DYNAMICS AT LOW 
DRIVE-TEMPERATURES (v < 0): 
SPATIO-TEMPORAL CHAOS 



We now investigate the low-temperature v < phase, 
in the absence of noise. It is convenient to work with 
dimensionless variables, obtained by rescaling x±, xu, t 



and S in Eq. (gj 



v\ V u 




B —> B 



uH d / 2 ~ 4 

Z 1/2 d/2-1 

or M r. 



leaving A as the only 'drive' parameter in the equation of 
motion, 



dS_ 
~dt 



dj + VIJS + S - (S ■ S)S - XS x d\\S 
gil Sx dfS + gxSx V 2 ± S . (32) 



We first investigate the static, spatially homogeneous 
steady states : 

(i) 'Paramagnetic steady state' represented by (S a ) = 
(average is taken over the steady state configurations) 
is a solution of the stationary equations. It is easy to see 
from Eq. © that this steady state is linearly unstable. 

(ii) 'Ferromagnetic steady state' with broken 0(3) 
symmetry represented by (Si) — (S2) — and (S3) = 1 
is also a solution of the stationary equations. This turns 
out to be linearly unstable too, as can be seen by per- 
turbing about this state by a small fluctuation u(x, t) (to 
avoid a clutter of terms we set g± = gu = since these 
terms are of higher order in gradients than the A term), 



d t ui(x,t) = V 2 mi(x, t) + A<9||U 2 (x, t) 
0tu 2 (x,£) = V 2 w 2 (x, t) - A<9||Ui(x, t) 
ftu 3 (x,t) = V 2 w 3 (x,<)-2w 3 (x,t). 



(33) 



Using the combination u + — u\ + iiL2, u = u\ — iu2 and 
113 the above equations simplify in Fourier space, 



dtu+(t) = -k 2 u+(t) + \k llU +(t) 7 
dtu u (t) = -k 2 u u (t)- AA||t£(t), 



-k 2 u 3 k(t) - 2u 3k (i) , 



(34) 



clearly showing that u^(t) = uj(0) exp(— k 2 t ± \k\\t) are 
unstable at large wavelengths when fcii < A. 

(iii) We next look for static, spatially inhomogeneous 
steady states; a natural candidate is the 'helical steady 
state' which is more conveniently expressed in terms of 
the variables p = \J S\ + S 2 and 4> = tan^S^/Si). In 
these variables, Eq. 10} for gn — g± — becomes 



dp 
at 

dip 
~dt 
dSs 

at 



V 2 P - p(V0) 2 +p-( P 2 + S 2 3 )p - XpSsd^ , 

V 2 + -(Vp) ■ (V0) + -(5 3 9||p - pd\\Ss) , 
p p 

V 2 S 3 + S 3 ~ (p 2 + S 2 )S 3 + x P 2 aa. (35) 
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A regular helix p = a, (f> = pxn and S3 — b (a, b and p 
are arbitrary constants) is a steady state solution if the 
projection of the local spins along the || axis b and the 
pitch 1/p satisfy the following relations 

2b 2 = 1 - a 2 (l + A 2 ) ± V(a 2 (A 2 + 1) - I) 2 - 4a 4 , (36) 
and 

2p = -Xb ± ^X 2 b 2 - 4(R 2 - 1) , (37) 



where R = yja 2 + b 2 is the magnitude of each spin. The 
only free parameter a is however bounded by a < (3 + 
A 2 ) -1 / 2 , coming from the requirement that b be real. 

Unfortunately even the helical steady state is linearly 
unstable as we show explicitly. Consider small fluc- 
tuations about the helical steady state, p = a + dp, 
4> = px\\ + c + 5(j) and S3 = b + u. To linear order the 
Fourier components of the fluctuations evolve as (for sim- 
plicity we exhibit the modes with fcj_ = 0) 




-k 2 - 2a 2 
ik\\(2p + \b) / a 
-2a(b - Xp) 



-ik\ia(2p 




(38) 




(C) 

FIG. 3: Plot of the real part of the eigenvalues ei, ei and ez 
versus a and fcii for A = 2. The eigenvalue e 1 is positive for 
all values of a and k\\ indicating linear instability. 



The signature of instability is that the real part of any 
one of the eigenvalues of the matrix D be positive. Fig. 
13 shows 2-dimensional plots of the real part of the eigen- 
values versus a and for a particular value of A. This 
shows that at least one eigenvalue has a positive real part 
for a continuous band of k\\ . We have checked that this 
result holds for other values of A. This implies that there 
is an infinity of unstable spatially periodic steady states 
parametrized by a (and for each value of a there are two 
values of b and p), a fact that will be of some significance 
later. 



The inhomogeneous helical steady state was suggested 
by the chiral nature of the driving. It would be a hard 
task to do an exhaustive check of all inhomogeneous con- 
figurations for possible steady states. Our strategy is 
therefore to solve the noiseless equations of motion (@J 
numerically starting from generic initial configurations 
in both 1 and 2 dimensions. The dynamics could either 
take the system to some other non-trivial inhomogeneous 
stationary state or lead to temporally periodic or chaotic 
configurations |27| . 

The numerical scheme for solving Eq. 0} should be 
chosen carefully as the linear derivative in the drive would 
give rise to numerical instabilities if the standard Eu- 
ler scheme of discretization were implemented |28j . We 
adopt an operator splitting method [23] which allows us 
to treat the dissipative terms and the drive separately 
under different discretization schemes. The dissipative 
part is solved using the standard Euler method [2£| and 
for the drive we use the following algorithm. The time 
evolution of the spins with the drive alone is a preces- 
sion about the local magnetic field /i(x, t) = d\\S(x, t). If 
/i(x, t) were a constant in space and time, the local spin 
S(xl, t) would have precessed about this field, keeping its 
magnitude |>S| fixed but changing its azimuthal angle cj) 
(taking the direction of h as the z— axis) by |/i(x, t)\At 
in a time interval of At. This would have been exact if 
the field h were a constant ; in our case however /i(x, t) 
depends on space and time and we introduce errors of 
0(A£). We choose At small enough so as to reduce this 
error. The advantage of this method is that it does not 
give rise to numerical instabilities and automatically pre- 
serves the magnitude of the local spin \S\ in time. In the 
simulation space and time are discretized with Air = 1 
and At = 0.0001 on a system of size N — 200 (large 
enough to avoid finite size effects) with periodic bound- 
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ary conditions. The local field is calculated by the rule 
hi = (Si+i — Si-i)/Ax. This field is used to update the 
local spin by the precession algorithm. 

Using this numerical scheme we can compute the time 
series of observables like the magnetization and energy 
E = J dx(V S) 2 . We first note that these quantities 
never seem to settle to a constant value, strongly sug- 
gesting that no stable steady state exists. The motion 
could therefore be either (quasi)periodic or chaotic. This 
should be revealed in a power spectrum analysis ; regular 
periodic motion would appear as a set of sharp delta 
functions. In Fig. 0^ we display the power spectrum 
P(oj) = \M3(uj)\ 2 of the third component of the total 
magnetization M3 = (S3) for data collected over more 
than 3 decades. The displayed spectrum has a clear 
smooth component with some small features which may 
be erased by more averaging and more sophisticated bin- 
ning. In addition, we find that the power spectrum fol- 
lows a power law (1/w 2 ) behavior over roughly 2 decades. 
The power spectrum of the total energy shows a similar 
behavior. These results strongly suggest that the dynam- 
ics is temporally chaotic 30]. We have also checked that 
this chaotic behavior shows up from a variety of initial 
configurations. The power spectrum we obtain from a 
study of our model in 2 space dimensions shows similar 
features (Fig.0j3). 

Since the components of spin obey partial differential 
equations (PDEs), we also check for spatial chaos. This 
is best visualized by constructing space-time plots of lo- 
cal quantities. For instance, Fig. [3^, is a space-time plot 
of the signed local pitch sgn(p) =sgn(9 x 0), suggesting 
spatio-temporal chaos |30j. 

We also evaluate, in d = 1, the space-time correlators 
of the spin C(x,t) = (5(x + x', t + t) ■ S(x ,<')) in this 
spatiotemporally chaotic phase. Fig. [5Jd shows that the 
equal-time spatial correlation function C(x, 0), decays ex- 
ponentially with a correlation length of the order of the 
lattice spacing. On the other hand, the unequal-time cor- 
relator C(0,t), seems to have a power- law over slightly 
more than a decade with an exponent 0.7. While this 
is not inconsistent with spatio-temporal chaos [30j . one 
generally expects power law correlations in non-critical 
systems only if a conservation law or a Goldstone mode 
is present. We do not understand the origin of this power 
law decay at present, 

Admittedly such characterizations of spatio-temporal 
chaos are only suggestive and should be made rigorous 
by studying the dependence of the number of positive 
Lyapounov exponents on system size. In spite of this, 
we hope we have provided convincing evidence that the 
asymptotic configurations in the low drive-temperature 
regime exhibit spatio-temporal chaos. The numerical ev- 
idence we presented was for d — 1 and 2, and though we 
cannot be sure whether this spatio-temporal chaos will 
persist at higher spatial dimensions, we feel that this is 
quite likely. This is because in our stability analysis of 
steady states done for arbitrary spatial dimension, we 
failed to find any reasonable stable steady state configu- 
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FIG. 4: (a) Log-log plot of y = \J uj 2 \Ms(lj)\ 2 versus oj in 
d = 1 showing the 1/uj 2 dependence of the power spectrum 
over approximately 1.5 decades, (b) Same plot for d = 2. 



ration at low drive-temperatures. Moreover, a Lyapunov 
stability analysis of the simpler equation dtS = \S x d\\S 
in arbitrary d reveals that a tiny disturbance in the initial 
conditions grows exponentially in time. Several questions 
arise, to which we do not have answers at present, such as 
whether there exists a low-dimensional chaotic attractor 
and if so what is its nature and dimensionality. 

The spatio-temporal chaotic phase that we just discov- 
ered has embedded in it an infinity of unstable (spatially) 
periodic steady states. It seems likely [jOlH^I that start- 
ing from generic initial conditions the configuration of 
spins would eventually visit these periodic steady states, 
although the time taken to visit any one of these peri- 
odic configurations is unpredictable. Since these periodic 
steady states are unstable, once visited, the dynamics will 
veer the spin configurations away from it. We now ask 
whether we can arrange that the spin configuration stays 
put in a prescribed periodic steady state having visited 
it ? This is the subject of control of spatiotemporally 
chaotic systems, one of the most important problems in 
modern chaos research [33] ■ There are two aspects to 
the control of chaos, stabilisation and targeting. Holding 
the periodic steady state having visited it is termed stabi- 
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FIG. 5: (a) Space-time plot of the signed local pitch, sgn(p) = 
sgn(dx<j}) (black and white patches denote sgn(p) = ±1 re- 
spectively), revealing spatio-temporal chaos. (b)Plot of equal- 
time correlation function C(x,0) vs x (see text). 



lization. However since the time taken for this visit from 
an arbitrary initial condition can be extremely large, it is 
desirable to target a prescribed unstable periodic steady 
state. There have been many proposals for controlling 
chaos in finite dimensional dynamical systems |3lL l32| . 
However there has been very little work in the more im- 
portant area of control of spatio-temporal chaos in PDEs 
(which correspond to an infinite dimensional dynamical 
system, see Ref. j33| for a review). Accordingly, it is sig- 
nificant that we are able to stabilize, target, and hence 
control spatiotemporal chaos in our model, as we now 
show. 



V. TARGETING AND CONTROL: 
EMERGENCE OF HELICAL STATES 

The helix solutions of Eq. I|35[l f° r v < are an infi- 
nite family of unstable spatially periodic steady states of 
the type discussed in [31| • These helical steady states are 
parametrized by a, the projection of the spin along the _L 
axis, p, the inverse pitch, and b, the projection of the spin 
along the || axis. Can chaos in our model be controlled 
so as to stabilize and target [31| these helical states? Our 




FIG. 6: Picture of the stabilized helical configuration from 
numerical simulations. 



control strategy focuses on the spin component 53. Thus 
for instance, in order to stabilize a specific helical config- 
uration (with fixed a, b and p), we could in principle wait 
till the dynamics (presumably ergodic) eventually leads 
to this configuration, after which we apply small pertur- 
bations to prevent S3 from deviating from the value b. 
This successfully stabilizes the prescribed helix, Fig. 

In order to target this prescribed helix, we add to i|35|) 
terms which would arise from a uniaxial spin anisotropy 
energy V 3 = r 3 (Sf-fc 2 ) 2 or V 3 = r 3 (S 3 ~b) 2 . We find that 
a sufficiently large and positive r 3 forces S 3 to take the 
value b exponentially fast starting from arbitrary initial 
configurations. The subsequent evolution, given by Eq. 
(I35|l on setting S3 = b is (in d dimensions), 



dp 
dt 



dt 



V 2 p ~ P (V0) 2 + p - (p 2 + b 2 )p - Xbpd\\ <t> , 



= V 2 0+-(Vp)-(V ( « + -(cV). 

p p 



(39) 



We now note that these equations can be recast as purely 
relaxational dynamics, 



dp 
dt 

at 



SF_ 
Jp-> 
1 SF 



(40) 



where the 'free-energy functional' F has the form of a 
chiral XY model |2Q|. 



/ - 2 I <!■'.,■ 



(Vp) 2 + p 2 (Vt) 2 ~(p 2 + b 2 ) 



-i( P 2 + b 2 ) 2 + Xbp 2 d ll 



(41) 



It is easy to see, using the chain-rule, that F is a Lya- 
punov functional of the dynamics |2l| , 



§-'■"■ 



d a x 



SF dSj. 
SS\~dT 

SF 



5F dS 2 
JS~2~dt 

SF 

s!h 



<0, (42) 



which decreases monotonically in time. Completing the 
squares, we see that du(p appears in F in the combination 
(l/2)p 2 (d\\cf> + \b/2) 2 , which is minimized by the helix 
4> = — (l/2)Abx||. Hence starting from any initial config- 
uration the system plummets towards the minimum of 
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FIG. 7: (a) Contour plot of the 'free-energy' Fh as a function 
of a and p. (b) Projection of plot (a) on the pa-plane. It is 
clear from the plots that there is a unique helical minimum. 



this F, which is a unique helix with parameters b, a and 
p (p is related to a and b via Eq. H37[H. That there is 
a unique minimum can also be seen by determining the 
'free-energy' Fh of the helical configurations from Eq. 
(|4l"|> and plotting Fh [p, a] against p and a (Fig. EJ . 

The stability of this 'free-energy' minimizing helix 
can also be tested directly from Eq. As before, 
we perturb about this helix : p = a + p(x, i), and 
= px\\ + c + 0(x, £), and deduce the growth of these 
perturbations to linear order, 



dp 
Of 
50 
0t 



V 2 p - 2a z p - a(2p + \b)d\\(j>, 
V 2 4>+-(2p + \b)d ]l p. 



(43) 



By going to Fourier space we may evaluate the eigen- 
values c^i(k) and uJ2(k) of the dynamical matrix, given 
by 



Wl>2 (k) = -(fcf! + fci+a 2 )± 



1/2 

k\{2p + Xb) 2 \ (.44) 



The helix configuration is stable if these two eigenvalues 
have negative real parts for all k, which they do, as can be 
seen numerically by substituting the values a and p take 
at the 'free-energy' minimum in the above expression. 

Let us now see whether our control is robust against 
noise. Indeed no discussion of the control of spatio- 
temporal chaos is complete without considering the ef- 
fects of noise which might result in occasional escapes 
from the otherwise well controlled system. Do these es- 
capes lead to an instability of the targeted configuration 
? We therefore modify l|43|l by including the nonconser- 
vative Gaussian white noises r\ p and 77^, to linear order. 
We shall declare the controlled helical state to be robust 
if the means (p) and (0) vanish and the variances are 
finite in the thermodynamic limit. We therefore ask for 
the statistics of small fluctuations with Fourier compo- 
nents pu(t) and 0k (t) about the controlled helical state, 
where 2ir/L < k < A for a system of linear extent L (A 
is the ultraviolet cutoff). It is clear from (|43|l that the 



means (p"k) and (0k) decay exponentially to zero: the re- 
laxation time for pk is finite at small fc, whereas that for 
0k goes as k~ 2 . 

To calculate the variances, note that the dynamics 
is governed in the mean by the Lyapounov functional 
(|41|l . and that the noise is spatiotemporally white. It 
follows [2(| that the steady-state configuration proba- 
bility P[p, 0] oc e~ cF where c is an effective inverse 
temperature. Since F ~ J[const(p) 2 + const(V0) 2 ] for 
small fluctuations about the helical minimum, i.e., P is 
approximately gaussian, this immediately implies that 
(|0k| 2 ) ~ k~ 2 and (|pk| 2 ) ~ const, for small k. Thus 
the variance (p 2 ) = J fc (|pk| 2 ) (2n/L < k < A), is L- 
independent for L — > 00 in any dimension d, whereas 
(^ 2 ) = Jfc(l^k| 2 ) diverges as L and InL respectively for 
d = 1 and 2, and is finite for d > 2. 

To obtain a more explicit asymptotic form for the the 
variance we calculate the equal-time correlation functions 
Cp(k,t) = {p k (t)p- k (t)) and C^(k,i) = (0 k (i)0_ k (i)) 
from the linearized equations of motion, 

^ = V 2 p - 2a 2 p - (2ap + \ab)d\\4> + r] p , 
<90 _,t 1 



V 2 0+-(2p+A6)5||p + ?/ , 



(45) 



where 77 p (x, t) — ?7i(x, t) cos(px||) + 772 (x, i) sin(pa;|| ) and 
T]<i,(x,t) = a _1 [77 2 (x, t) cos(px\\) - 771 (x, i) sin(pcc|| )]. The 
noises rj p {k, t) and 770 (k, t) satisfy 



(77 p (k,i)77 p (k',<')) = 2fl<5 k ,_ k ,<y(t - t') . 



2B 



6 K -u>6(t - t') . 



(^(M^k'.O) = 0. 



(46) 



We may easily solve Eq. (|45fl for Fourier transforms pk 
and 0k, 



0k M 



n,, 1 (i] p {k,uj)(iuj + k 2 
D{k,uj) 

-irj^ik, uj)k\\a(2p + Xb)) . 
1 



D(k,LU 



■(rj^Lujiiuj + k 2 + 2a 2 ) 



-77 p (k,w)fc||(2p+A6)), 



(47) 



where 



D(k,uj) = -uj 2 +2iuj(k 2 +a 2 )+k 2 (k 2 +2a 2 )-k' 2 (2p+\bj 2 . 

(48) 

We now compute the equal time correlation functions 
Cp and Ct averaged over the noise: 



(49) 



C p (k) = (pk(t)p-kW) = J (Jw(pk(w)p-k(-w)) 

k 4 + k 2 a 2 (2p + Xb) 2 B 
= B 2/+/_(/+ + /_) + 2(/ + + /_) 
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where /+ and /_ are the roots of the equation D(k, ui) 
0, 



a 4 + k 2 (2p + Xbf 



-l 1/2 



(50) 



and 



c^(k) = (Mm-k(t)) = J $cj(M")<l>-k(-u)) 

B k 2 + 2a 2 + k?,(2p + Xb) 2 



B 



2/ + /-(/ + +/-) 
1 



2a 2 /+ + f- 



(51) 



We use the k — > behavior to evaluate the variances 
Ap and At. In this infrared limit 



C p (k,t) 
Ci(k,t) 



k 2 



B 



kl 



k 2 \ 1 



(52) 



The variances Ap and At obtained on integrating the 
correlators over all k and then taking the thermodynamic 
limit L — > oo depend sensitively on the spatial dimension 
d. The variance of p is finite in all dimensions, 



7r /a d = 1 

ln(l + (A/a) 2 ) d=2 
finite d = 3 



(53) 



while the variance of <f> diverges in 1 and 2 dimensions 
and is finite in higher dimensions, 



L d=l 
In L d=2 
finite d = 3 



(54) 



Thus occasional excursions from the controlled state as 
a result of the noise do not lead to an instability of the 
targeted state for d > 2; the behavior for d < 2 is no 
worse than for a thermal equilibrium XY model. 



VI. DISCUSSIONS AND FUTURE DIRECTIONS 

In this paper we have constructed the simplest exam- 
ple of a spatially extended dynamics in which dissipation, 
(reversible) precession and spatially anisotropic driving 
act in concert to produce spatio-temporal chaos in a 
whole region of parameter space. The model, the classi- 
cal Heisenberg magnet in d space dimensions, is driven 
by imposing a background steady current of heat or par- 
ticles (or any other mobile species) which couples to the 
Heisenberg spins. Our driven Heisenberg model (DHM) 
is a natural generalization of the DDLG models to the 



case of a 3-component axial-vector order parameter S, 
and is an important step in the exploration of dynamic 
universality classes f ar from equilibrium 0, E3 ■ 

In the limit where the mobile species is 'fast', the im- 
posed steady current alters the effective dynamics of the 
isotropic magnet so as to generate a driven precession 
term XS x d\\S, that is responsible for all the remarkable 
phenomena we predict. We summarize the nature of the 
asymptotic configurations in a 'non-equilibrium phase di- 
agram', Fig. |H1 as a function of the drive-temperature 
v. The system exhibits a 'paramagnetic steady state' at 
high v and a 'critical steady state' at v = Q, with power- 
law correlations induced by the driving even when the 
equilibrium Heisenberg magnet (without the driving) is 
paramagnetic. The drive takes the system away from the 
Wilson-Fisher fixed point leading to a new drive induced 
universality class. At low drive-temperatures both the 
homogeneous and inhomogeneous steady states are un- 
stable. In particular, the system has an infinity of spa- 
tially periodic unstable steady states which are helical. 
We have provided evidence, at least in d = 1 and 2, that 
the dynamics at low v is spatio-temporally chaotic. We 
have found that the spatio-temporal chaos may be 'con- 
trolled' to target any desired helical steady state. This 
control works even in the presence of noise in dimensions 
d > 2. The generic occurence of spatio-temporal chaos 
in our model is encouraging, given the intriguing formal 
similarity to the nonlinearity (v- V)v in the Navier Stokes 
equation. 

Our criterion for declaring the low temperature phase 
as spatio-temporally chaotic is based on rather well estab- 
lished tests of power-spectrum analysis, space-time plots 
and decay of space-time correlators. Admittedly such 
characterization is not very rigorous; ideally one would 
like to analize the system size L dependence of the Lya- 
pounov spectrum. We leave this for future investigation, 
as also several questions regarding this spatio-temporal 
chaotic regime, such as the existence and nature of a 
low dimensional attractor or the possibility of complex 
ordered states in higher dimensions, 

Note that while the entire low drive-temperature phase 
is chaotic, the critical phase is not; the DRG analysis def- 
initely shows a stable fixed point in all dimensions. This 
is quite different from the possibility of a chaotic crit- 
ical point, discussed in Ref. Q, for a general class of 
anisotropically driven 0(n) models along du directions. 
Their driving appears as a noncquilibrium source of noise 
which is not bound by a fluctuation-dissipation relation. 
At the critical point of this model they find, to 1-loop, 
no stable fixed point when d\\ > df,. They tentatively 
put forward the possibility of spatio-temporal chaos at 
criticality. There have been earlier suggestions of chaotic 
criticality in equilibrium critical behavior, albeit in sys- 
tems with quenched disorder [33| . 

We conclude this paper by reminding the reader again 
of the variety of explicit lattice and continuum realiza- 
tions of our driven dynamics discussed in Sect.IIB . We 
hope that this will stimulate a search for experimental 
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FIG. 8: Non-equilibrium phase diagram as a function of the 
drive-temperature v. 



systems, e.g., isotropic magnets carrying a steady parti- 
cle or heat current, as well as model magnetized lattice- 
gas simulations, where the predictions of our model can 
be tested. 
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APPENDIX A 

In this Appendix we present the details of pcrturbative 
calculation mentioned in Sect. IV. The diagrams corre- 
sponding to the lowest order terms in the perturbation 
expansion are constructed following the rules shown in 
Fig. El 

Corrections to Go(k, u) 

(I) Corrections from the A vertex : 

Graphs (I) and (II) in Fig. 1101 show the one-loop cor- 
rections to Go(k, ui) due to the A vertex. 



SA(k,w) 



2X 2 B 



4 



(<?|| - (fe|| - <?||))( fc || - (-<?||))G (k- q,w - is)C (q,v) 



+((k\\ - q\\) - q\\)(k\\ + (fc|| - q\\))Go(q,u)Co(k -q,w-v) 
(2g||-fcy)(*5|| +«||) (fc||-2g||)(2g||- 9 ||) 



7 2 (q) 



7 2 (k-q) 



7 2 (q)+ 7 2 (k-q) 



(Al) 



r 



In the above integrals 7 ™(k) is defined as 



tribution coming from the A vertex to order k° as (£q)a> 



7 ™(k) = 



r\\kf\ + rik 2 , 



(So); 



X 2 B 



fiq 



7 4 (q) 



ttAX 2 B 

' , 3/2 1/2 

4r,i r 



(A2) 



We expand £,\(k, 0) in powers of fc|| and k±. Terms which 
are higher order than 0(k 2 ) or 0(k 2 _) are ultraviolet con- 
vergent. Coefficients of the terms of order k°, k 2 and k\ 
denote changes to v, rii and r± respectively. 

(i) Term proportional to k° 

This part of Sx(k, 0) does not have any infrared diver- 
gences but has divergences in the ultraviolet. We intro- 
duce an upper momentum cutoff A to get an estimate of 
the correction to v. We will denote the self energy con- 



(ii) Term proportional to k 2 



( S ll) 



where 



X 2 Bk 2 



dq 



X 2 Bv e / 2 r ]l kf l 
8n 2 e(r ll r ± ) 3 / 2 ' 



ctq 



16r 



f/2 



7 4 (q) 87r 2 (r||r ± ) 3 / 2 



2 7 4(q) 7 8( q ) 7 6 (q ) 



(A3) 



(A4) 
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a a 



G (k, a) ) 



k , to 



q ,v 



k - q , co - v 



7 q 2 .v 2 



C (k,co ) 



2 W < k if i„M 



= - uF. 



where we have evaluated the integral 



(t q 



v -e/2 



q \\ q 



7 8 (q) 487r 2 rx(r||rj_) 3 / 2 

q , v 




(A9) 



k , to 



k - q- q 2 , m-V,- V 2 



FIG. 9: Graphs for correlation functions and vertices at the 
tree level. F a p~ /S = (l/3)[£ aj g<5 7 4 + 8^5^ + 5 a */Sps]- 



q,V 



k- q , (o-v 





— k,co \ / — k,co — k,co \ J — k,(B 

k - q , to- v q i v 

( I ) ( II ) 

FIG. 10: 0(A 2 ) corrections to Go(k,w) 



tl (7 



e/2 



1 



7 6 (q) 327T 2 (r||r ± ) 3 /2 



7 8 (q) 64 7 r 2 (rfi /2 r5 /2 ) 



(A5) 



(A6) 



have been evaluated using the standard integrals for the 
general form J dq(q^q b j_)/j c (q) given in Ref. |2^. We 
have also had to make use of the asymptotic expansion 
for the Gamma function r(— n + e), when n is zero or any 
positive integer and e — > 0, 



r(-n + e) 



(-1)" 



1 

1+ 2 



1 

7 

n 



where 7 is the Euler-Mascheroni constant [3J 



0(e) 
(A7) 



(iii) Term proportional to k 2 



(Sx)a = - 



X 2 Bk 2 , 



qq \\ 



16rj_c7j 



4rj_ 



7 8 (q) 7 6 (q) 



5X 2 Br ± k 2 ± v- e / 2 1 
967r 2 rj_(r||r'j_) 3 / 2 e 



(A8) 



k , co 



k , co 



FIG. 11: 0(u 2 ) corrections to the k = part of Go(k, 



k - q 



, co-v 




k , co 



FIG. 12: 0(A 2 ) corrections to C (k,u). 



(II) Corrections from the u— vertex : 
There is also a correction to the propagator G(k, oS) 
coming from the u— vertex; to 0(e) this is shown in the 
figure below. 

To 1-loop, the correction to G(k, u>) from this interac- 
tion comes only from the k = piece given by, 



(E„)t 



5u 

y 

5itAuB 



S qS ^Go(q, v) 
1 



1/2 3/2 



(A10) 



The net 1-loop correction to the propagator G adds up 



to 



E(k,w) = [(S ) A + (E )„] + (S||) A fc 2 + (S ± ) A fci. (All) 
Corrections to Go(k, w) 
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"T / s Q Sl/ (( k \\ - - - ( k \\ - 9||))Co(k- q,w - v)C {(\,v) 



J 7 6 (q) 

~ 327r 2 (r||r_L) 3 /2 e ' 



(A12) 
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FIG. 13: 0(A 3 ) correction to the A— vertex. 



The renormalized couplings are obtained from the ver- 
tex corrections to A and u. 



Corrections to the A— vertex 
The corrections to A come from the three-point correla- 
tors, 



The contribution from (I) to the vertex function T\ is 
given by, 



-i\ 
~2~ 



((fci|| ~ <l\\) + ki\\){{k 2 \\ +q\\) + 9||)((fc 1 || +fc 2 ||) + (fc 1 || -q\\)) 



Go(q,^)G (k 2 + q,w 2 - f)Co( k i ~q> w i - ^) 

a 3 s /■ (8fcin + *2||)<jf 



16 



7 6 (q) 



(A13) 



Likewise the contribution from (II) to T\ is 



-iX 
~~2~ 



ct qdv 



(-fci|| -9||)((*2|| +9||) + «||)((fei|| +A2||) + (*i|| -«||)) 



G (ki - q, wi - ^)G (k 2 + q, w 2 - i/)C (q, f) 

iA 3 S /■ (2*i|| +^ 2 ||)g 2 
ag- 



32 



7 6 (q) 



(A14) 



while (III) gives 
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(-(fan +0||)-(&i|| +*2||))((*2|| +9||)+9||)((*i|| -9||) + Ai||) 

G (q, i/)G (ki - q, wi + f )G (k 2 + q, Wi + v) 
(2fc!|| -5fc 2 ||)«f 



»A 3 g 
32 



cf g- 



7 6 (q) 



(A15) 
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FIG. 14: 0(tt 2 ) correction to the u— vertex. 

The three contributions combine to give 

3iX 3 B ( 9y(2fci|| + k 2 



r A (2fcin +fc 2 ||) 



16 7 7 6 (q) 

3iA 3 B(2fc 1 || +fc 2 ||)v- e/2 
5127r 2 (r||r ± ) 3 / 2 



|A16) 



Corrections to the u— vertex 
There are two contributions to the vertex function T u to 
0{u 2 ) and 0(A 4 ). 

The 0(u 2 ) graph gives 



= F, 



Hit 2 f 

F af 3 7 s—g- / dqdvC (q, i/)G (ki + k 2 - q, Wi + w 2 - i/) 

(A17) 



llu 2 S 

11m 2 B 1 



The 0(A 4 ) graph gives the following contribution : 



J 



-1152F a p lS (J^j j SqSv{{k n - - g||)(fc 4 || - {k n - - q\\))(k 3 \\ - {k 2 \\ + q\\))(k 2 \\ + q\\) x 
G (ki q, uj - i/)G (ki - k 4 - q,wi - w 2 - f)G (k 2 + q, w 2 + i/)C (q, 



= — F, 



27A 4 



aft-yS 



(A18) 



The net correction to the w— vertex is, 

11m 2 £ 1 27A 4 



ll7r 4 K 2 27tt 4 t 2 



12Se 

where k and r are defined by 

uB 

K ~ 2TrHr\ /2 rT) 



8Be 



(A19) 



(A20) 



XB 



T = 



2^{rl /2 r 3 / 2 ) 
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